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Introduction 

In  this  paper,  we  give  a  theory  of  spectral  approximation 
for  closed  operators  in  Banach  spaces.  We  then  apply  this  t:  eory  to 
an  analysis  of  the  approximation  of  the  spectral  oroperties  of 

some  differential  systems  by  finite  element  methods. 

Bramble  and  Osborn  [1]  and  Osborn  [14]  developed  a  theory  of 
spectral  approximation  for  compact  operators  in  Banach  spaces. 

Their  theory  can  be  applied  to  the  analysis  of  many  numerical 
procedures  for  the  spectral  approximation  of  differential  operators, 
T  ,  such  that  T  +  XI  has  a  compact  inverse  for  some  X  6  I  . 

Most  of  the  differential  systems  in  the  theory  of  elasticity  are 
in  this  class. 

However,  there  are  many  differential  systems  of  interest  in 
mathematical  physics  which  do  not  have  compact  resolvents.  These 
operators  can  have  continuous  spectrum,  eigenvalues  of  infinite 
multiplicity,  and  finite  limit  points  of  eigenvalues.  Also,  the 
eigenfunctions  need  not  be  smooth  since  the  differential  systems 
are  not  necessarily  elliptic. 

Descloux,  Nassif,  and  Rappaz  [4,5]  have  studied  the  approxi¬ 
mation  of  the  spectrum  of  a  differential  system  of  interest  in 
magnetohydrodynamics  which  has  a  bounded  inverse,  but  not  a 
compact  inverse.  They  developed  a  theory  of  spectral  approximation 
for  bounded  operators  which  treats  this  problem.  An  analysis  of 
the  approximation  of  the  spectral  properties  of  a  class  of  bounded 
operators  by  finite  element  methods  has  also  been  done  by  Mills 
[12,  13]. 
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The  results  in  this  paper  apply  to  closed  {not 
necessarily  bounded)  operators  in  Banach  spaces.  We  apply  the 
perturbation  theory  developed  in  this  paper  to  the  study  of  a 
finite  element  procedure  for  approximating  the  spectral  properties 
of  a  differential  system  modeling  a  fluid  in  a  rotating  basin. 

We  note  that  unlike  previous  authors ,  we  analyze  the  approxi¬ 
mation  of  the  differential  operator  directly  and  not  through  its 
inverse. 

The  time  dependent  equations  for  the  differential  system 
modeling  a  fluid  in  a  rotating  basin  are 


(x,t)  e  a  x  ®  , 


(0.1)  ||  *  -7?  -  f$  -  wRu  ,  ( x ,  t )  e  (2  x  TR  , 

u  •  n  =  0  ,  (x,t)  e  3!1  *  1  , 

/  £dx  =  0  , 

$2 

2 

where  (2  C  ®  is  a  bounded,  connected  open  set  with  smooth 
boundary,  3(2  ,  u  *  (u^,U2)  is  the  horizontal  volume  transport, 

Z  is  height  of  the  fluid  above  equilibrium  level,  R  is  the 
linear  operator  Ru  ■  (-U2,u^)  ,  n  is  the  exterior  normal  to 
3(2  ,  and  f  _>  0  and  oi  are  real  constants  representing 
friction  and  Coriolis  terms. 

Thus,  we  are  concerned  with  approximating  the  spectral 
properties  of  the  system 

T(C,u)  =  (-7  •  u,  -  7?  -  fu  -  u)Ru)  ,  x  e  (2  , 

(0.2)  u  •  n  =  0  ,  X  e  3(2  , 

/  cdx  **  0  . 

(2 

We  note  that  T  -  I  is  a  symmetric,  formally  dissipative  operator 
with  maximal  non-positive  boundary  conditions  [9]  . 
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If  f  =  a)  =  0  ,  then  T(0,RV>l>)  =  0  for  all  smooth  functions 
ip  such  that  4>  *  0  on  .  Hence,  in  this  case  (f  =  cj  =  0)  , 

0  is  an  eigenvalue  of  T  of  infinite  multiplicity.  Since 
the  sum  of  an  operator  with  non-compact  inverse  and  a  bounded 
operator  has  a  non-compact  inverse,  it  follows  that  T  -  XI 
has  a  non-compact  inverse  for  the  general  case  f  ,uj  e  JR  ,  X  e  (E  , 

X  belonging  to  the  resolvent  set  of  T  . 

We  give  here  error  estimates  for  a  finite  element  procedure 
proposed  by  Platzman  [15]  to  approximate  the  spectral  properties 
of  (0.2)  .  The  self ad joint  case  f  *  0  was  analyzed  by  Luskin 
[11]  by  techniques  different  from  those  used  here. 

In  Section  1  ,  we  give  a  general  theory  for  the  approximation 
of  a  closed  operator  A  by  a  family  of  finite  dimensional  operators 
{A^}  .  In  our  applications,  A  will  be  a  differential  operator 
and  A^  will  be  an  approximation  of  A  given  by  a  finite 
element  procedure.  We  propose  two  properties  PI)  and  P2)  and 
show  that  these  properties  imply  the  convergence  of  the  spectral 
properties  of  A^  to  those  of  A  .  Error  estimates  in  the  applica¬ 
tions  will  follow  from  Theorem  1.3. 

We  give  special  results  in  Section  2  for  the  case  A  and 
A^  are  self ad joint  operators  in  a  Hilbert  space.  These  results 
apply  to  the  approximation  of  continuous  spectrum. 

The  operator  theory  developed  in  Sections  1  and  2  is 
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spectral  properties  of  differential  operators  by  finite  element 
methods.  In  Section  3  ,  we  apply  our  theory  to  obtain  error 
estimates  for  the  approximation  of  the  spectral  properties  of 
scalar,  second  order,  uniformly  strongly  elliptic  operators 
by  the  standard  finite  element  method.  This  example  is  included 
even  though  the  results  are  not  new  since  we  believe  that  its 
inclusion  will  make  it  easier  for  the  reader  to  understand  our 
main  example  in  Section  4. 


In  Section  4,  we  define  and  analyze  the  approximation  of  T  . 

Theorem  1.3  gives  optimal  order  estimates  for  the  approximation 
of  eigenspaces.  Optimal  order  eigenvalue  estimates  for  this 
problem  have  been  given  in  [11]  for  the  self-adjoint  case  and 
are  derived  from  the  results  in  this  paper  and  [3]  for  the 
general  case. 

We  note  that  we  use  without  explicit  reference  the  classical 
spectral  theory;  see  for  example  Kato  [6],  Riesz-Nagy  [16]  .  ^ 

i 


i 

I 

i 
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1.  Approximation  of  the  Spectrum  of  Closed  Operators  in  Banach  Spaces. 

We  first  introduce  some  notation.  Let  X  be  a  complex 
Banach  space  with  norm  ||  ||  .  Denote  by  (&  (X)  the  set  of 
bounded,  linear  operators  B:  X  •»  X  .  Also,  denote  by  C(X) 
the  set  of  closed,  linear  operators  C:  £^(C)  CX  +  X  where  the 
domain  of  C  ,  & (C)  ,  is  not  necessarily  dense  in  X  .  For 

C  e  C{X)  ,  p (C)  is  the  resolvent  set  of  C  defined  by 

pKc)  =  (z  e  <c  !  (z-c)"1  e  &(x)}  . 

If  z  e  P (C)  ,  we  define  the  resolvent  operator  R  (C)  *  (z-C)  : 

z 

X  -*■  X  .  The  complement  of  p  (C )  is  a  (C)  =  { z  e  <E |  z  0  p  (C) }  , 
the  spectrum  of  C  . 

Let  Y  and  Z  be  closed  subspaces  of  X  and  x  e  X  . 

We  set 

<5<x,Y)  »  inf  ||  x-y  ||  ,  <5(Y,Z)  *  sup  <5(y,Z)  , 

y€Y  yeY 

ifyll-i 

6(Y,Z)  =  max(<5  (Y,Z)  ,6  (Z,Y)  J  ; 

A 

6  (Y, Z)  is  called  the  gap  between  Y  and  Z  and  is  a  measure  of 
the  "distance"  between  these  spaces.  If  C  and  D  are  in  ^  (X) 
with  graphs  GC,GD  CX  *  X  ,  then  we  define  5(C,D)  ■  5(GC,GD)  ,  i.e., 

6  (C ,D)  »  sup  inf  {  ||  x-y ||  +  ||Cx-Cy||  }  . 

xe^c)  ye-d'(D) 

11*11+  ||CX  ||=  1 

Furthermore,  we  define 

5(C,D)  »  max(5  (C,D)  ,5  (D,C))  . 
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Finally,  if  X  is  subspace  of  & (C)  H  (d)  ,  we  set 

II  C-D|(  =  sup  ||  Cy-Dy (I  . 

yeY 

l|y|l-i 

Now,  let  A  e  ^(X)  be  a  given  operator.  In  order  to 
approximate  a (A)  ,  we  consider  a  family  {X^}  of  finite 
dimensional  subspaces  of  X  parametrized  by  h  and  linear  operators 

Ah:  Xh  -  Xh  ' 

We  denote  by  a  (A^)  ,  p  (A^)  ,  and  Rz  (Ah>  :  Xh  -*•  X^  the 
spectrum,  resolvent  set,  and  resolvent  operator  of  A^  considered 
as  a  bounded  operator  in  X.  .  However,  when  used  in  connection 

il 

with  expressions  of  the  type  ”&(A^,A)"  or  "  U  A— A^|  |  x  "  ,  Ah  is 

h 

considered  as  a  closed  operator  in  X  with  non  dense  domain  Xh  . 

Let  r  C  p (A)  be  a  given  Jordan  closed  curve;  then 

E  *  f  Rz  (A)  dz:  X  *  X 

and  if  r  c  P  (Ah)  , 

Eh  =  WI  /,  Rz(Ah)dz:  Xh  Xh 

are  the  spectral  projectors  relative  to  A  and  A^  respectively. 

We  shall  also  use  the  relations 

AE  *  —r-  /  zR  (A)  dz :  X  -  X 

2tt  l  z 

and 

Vh  *  sr  {  xh  *  xh  • 
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We  now  introduce  the  following  two  properties: 


PI) 


lim 

lr"0 


<5  (A^A) 


0  , 


P2)  Vx  e  X  ,  lim  <5  (x,X.  )  =  0  . 

h-0  h 

The  following  theorems  contain  the  main  results  of  this  section. 
We  note  that  c  <  ®  shall  denote  in  this  paper  a  positive 
constant  which  is  independent  of  h  ,  but  which  varies  from 
estimate  to  estimate. 


Theorem  1.1. 

Suppose 

PI)  is 

valid  an! 

let 

K  C  p  (A) 

be  a 

compact  set. 

Then  there 

exists 

h0  >  0 

and 

c  =  c(K) 

such  that 

for  h  <  hQ 

we  have  K 

C  p (Ah) 

and 

||Rz(A)-Rz(Ah)||x  <  c6(Ah,A)  ,  zeK. 

h 


Theorem  1.1  shows  that  if  PI)  is  valid,  then  the  approxima¬ 
tion  of  o  (A)  by  ct  (Aj^)  is  upper  semicontinuous .  Furthermore, 

for  h  <  hQ  ,  1 1  Rz  <Aj^)  1 1  x  is  uniformly  bounded  on  K  .  This 

h 

is  a  stability  property. 

Theorem  1.2.  Suppose  PI)  is  valid.  Then  there  exists  h^  >  0 

and  c  so  that  for  h  <  tig  we  have  the  bounds 

(1.1)  !!  E-Eh||  x  +  ||AE-AhEh||x  <c5(Ah,A)  , 

h  h 

(1.2)  6  (Eh  (Xh)  ,E(X)}  <  c5  (Ah,A)  , 

6  (x,Eh(Xh))  £  c{<5  (x,Xh)  +  5  (Ah,A)  ||  x||  }  ,  xeE(x)  . 


(1.3) 
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Theorem  1.3.  Suppose  that  PI)  and  P2)  are  satisfied  and  that 
E (X)  is  the  finite  dimensional  subspace  corresponding  to  an 
isolated  eigenvalue  X  of  algebraic  multiplicity  m  of  A  . 

Let  a  be  the  ascent  of  (X-A)  and  let  f  be  a  holomorphic 
function  defined  in  the  neighborhood  of  X  .  We  set 

Yh  =  mini* <A|E(X) ,Ah).S(Ah,A)f 

Then  for  h  small  enough,  A,  I  .  :  E,  (X,  )  *  E,  (X.  )  has 

h'E.  (X.  )  h  h  h  h 

h  n 

exactly  m  eigenvalues  X^  ^,...,Xm  ^  repeated  according  to 
multiplicity.  Also,  there  exists  hQ  >  0  and  c  such  that  the 
following  bounds  are  valid  for  h  <  hQ  : 


(1.4) 

6 [Eh(Xh) ,E (X) )  £  cyh  , 

(1.5) 

^AtJEh(Xh)  ,AIe(X))  -CYh  ' 

(1.6) 

l£(i>  '  5  X  £<h,h>l  i  ' 

(1.7) 

lx-xi.hl°ielrh  • 

Remarks.  If  PI)  and  P 2 )  are  satisfied,  then  (1.3)  of 

Theorem  1.2  shows  that  lim  5fx,E,  (X,))  =  0  ,  Vx  e  E(X)  ,  i.e., 

h-0  h  h 

any  x  e  E (X)  can  be  approximated  by  vectors  in  E^CX^)  .  It 
is  not  true  however  that  any  z  e  a  (A)  can  be  approximated  by 
eigenvalues  of  A^  ,  i.e.,  there  may  exist  z  e  a (A)  so  that 
dist(z,a (A^) )  fails  to  converge  to  zero.  (the  approximation 
of  a (A)  by  a (A^)  is  not  necessarily  lower  semicontinuous) . 
The  classical  counter-example  is  the  shift  operator  (see,  for 
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example,  Kato  [6,  p.  210]).  Note  that  if  A  and  A^  are 
selfadjoint  operators  in  a  Hilbert  space,  then  lower  semicontinuity 
is  true  under  weaker  conditions  than  PI),  P2)  (see,  for  example, 
Kato  [6,  p.  431])  . 

In  order  to  prove  Theorems  1.1  and  1.2,  we  first  prove  a 
sequence  of  lemmas  which  are  variants  of  results  found  in 
Kato  [6,  p.  197-208]. 

Lemma  1.1.  Let  B  8$(X),  C  e  <f(X)  .  Then 

a)  6  (C ,  B)  <  ||C-Bi|  £,(c)  , 

b)  IlC-BlI^cj 

<  (1+  ||  B  ||  )2<S(C,B)/{1  -  (1+  ||  B||  )  5  (C  ,B)  ) 
if  the  denominator  is  positive. 

Proof.  Part  a)  follows  directly  from  the  definitions.  We  prove 
part  b)  .  Let  x  e  fr'iC)  ,  ||x||  =  1  and  let  e  >  0  be 
arbitrary.  It  follows  from  the  definition  of  5 (C,B)  that 

there  exists  y  €  X  such  that 

(1.8)  (|  x— y  1 1  +  ||  Cx-By  ||  <_  5  (C,  B)  {  ||  x||  +  ||Cxj|}  +  e. 

Consequently,  we  have 

II  (C-B)x||  £  ||Cx-By||  +  ||  By-Bx|| 

<  ||Cx  -  By  1  (  +  ||  B 1 1  ||  y-x|| 

(1*9)  i  (1+  ll  B 1 1  ){  ||  Cx-By ||  +  II  y-x||  } 

(1+  ||B||  )<5  (C,B){  1 1  x  1 1  +  1 1  Cx  1 1  }  +  e  (1+  ||B||  )  . 


< 
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Replacing  (|  Cxjj  by  jj  (C-B)xj[  +  ||  B ||  ||X||  in  the  right 

hand  side  of  the  above  inequality  and  letting  c  ■*  0  yields 

the  estimate 

||  (C-B) x||  < 

(1.10) 

(1+  ||  B||  )$  (C,B){  (1+  ||  B||  )  ||  x||  +  |!  (C-B)xj|  }  . 

The  result  follows  directly  from  (1.10)  Q.E.D. 

Remark .  It  follows  from  Lemma  1.1  that  if  A  e  $  (x)  ,  then 

property  PI)  is  equivalent  to  lim||A-Av||  v  =  0  .  This  is  the 

h-0  -1  *h 

spectral  approximation  condition  of  Descloux-Nassif-Rappaz  [4]  . 

Lemma  1.2.  Let  B  6  <3  (X)  ,,.*d  C,D  e  £  (X)  .  Then 

6  (C+B,D+B)  <  (1+  ||  B|j  )26(C,D)  . 

Proof.  Let  x  e  ^(C)  ,  ||  x||  +  ||  (C+B)x||  =  1  and 

let  e  >  0  be  arbitrary.  We  have  by  the  triangle  inequality  that 

(1.11)  |f  X  (I  +  ||CX||<||X||  +  ||  (C+B)  x||  +  ||  Bx  ||  <  1  +  ||  B  ||  . 

It  follows  from  the  definition  of  6  (C,D)  that  we  may  choose 
ye  & (D)  such  that 

(1.12)  |(  x-y  (|  +  ||  Cx-Dyi|  <  5  (C ,D)  (1+  ||  B  ||  )  +e  . 

Hence, 

||  x-y  ||  +  ||  (G+B)  x  -  (D+B)y|| 

<  1 1  x-y |j  +  ||  Cx-Dy  [|  +  ||  B||  ||x-y|| 

(1+  ||  B ||  )  2S  (C,D)  +  (1+  ||  B ||  )  e  . 


(1.13) 


< 
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The  result  follows  directly  from  (1.13)  after  letting  e  -*•  0  . 

Q.E.D. 

Lemma  1.3.  If  C  and  D  S  £(X)  are  invertible,  then 

5(C,D)  =  MC-1,©"1)  . 

Proof .  This  result  follows  directly  from  the  definitions.  Q.E.D. 

Lemma  1.4.  Let  C,D  e  *£(X)  and  suppose  there  exists  <  such 
that 

(1.14)  ||  Dx  ||  >  <  ||x||  x  6  S'  (D)  . 

If  <5 (C ,D)  <  min(l,<),  then  C  is  invertible. 

Proof.  We  show  that  if  C  is  not  invertible,  then  <S(c,D)  _>  min(l,c 
Let  x  P  ^(C)  be  such  that  ||x||  *  1  and  Cx  =  0  .  Let 
e  >  0  be  arbitrary.  It  follows  from  the  definition  of  6(C,D) 
that  since  Cx  =  0  we  can  choose  y  e  (D)  such  that 

||x-y||  +  ||  Dy  ||  <  5  (C,D)  +  £  . 

It  follows  from  (1.14)  that 

6  (C,D)  +  e  >  ||x-y||  +  <  Jl y 1 1 

l  |i-  II y H  I  +  <  llyll 

^1,15)  _>  min  (1,<  )  { 1 1-  ||  y  ||  |  +  ||y||  } 

>_  min(l,<)  . 

Since  £  >  0  was  arbitrary,  it  follows  that  we  have  reached 
the  contradiction  <5(C,D)  >  mind,*)  .  Q.E.D. 
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Proof  of  Theorem  1.1.  In  this  proof,  c  depends  on  K  but  all 
estimates  are  uniform  for  z  e  K  .  By  Lemma  1.2, 

(1.16)  <5  (z-A.  ,z-A)  <  C(S  (A,  ,A)  ,  z  6  K  . 

n  —  h 

It  follows  from  (1.16)  and  Pi)  that 

(1.17)  lim  6  (z-A.  ,  z-A)  -►  0  uniformly  for  z  e  K  . 

h-*-0  h 

Since  K  c  p (A)  ,  there  exists  c^ ■  such  that 

(1.18)  ||  (z-A)x||  >_  ^  ||x||  ,  x  e  ^(A)  ,  z  e  K  . 

We  can  now  conclude  from  (1.17),  (1.18)  and  Lemma  1.4  that  there 
exists  tig  >  0  such  that  z  -  A^  is  invertible  for  h  <  hg  and 
z  6  K  .  Since  X^  is  finite  dimensional,  we  have  that  K  C  p (Ah) 
for  h  <  hg  . 

Furthermore,  it  follows  from  (1.16)  and  Lemma  1.3  that 

(1.19)  6(Rz(Ah) ,Rz(A))  <  c5(Ah,A),  z  6  K  . 

Hence,  we  can  obtain  from  Lemma  1.1b  the  result 


(1.20) 


II  Rz  (a>-r2  (Ah>  H  Xh  1  c5(R2(Ah).-Rz(A)J 

_<  c6  (Ah,A)j  h  <  hg,  z  e  K  . 

Q -E  .  D . 


Proof  of  Theorem  1.2.  The  result  (1.1)  follows  from  Theorem  1.1 


and  the  following  estimates: 
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(1.21) 


HE“EhlU  +  ll^'Vhll 

n  n 


-  list  ^Rz(A,-Rz(Ah»dsll  x„ 
r  n 

*  ii 2#r  /  i2Rz(A)  -  2  Rz(Ah))<i2ii  x 
r  h 

<  <  2tt  )  “ 1  max  ||  R  (A)  -R  (A.  )  ||  y  /  (l+[z|)|dz 

zer  2  2  n  r 


The  estimate  (1.2)  now  follows  directly  from  the  estimate  (1.1). 

In  order  to  prove  (1.3),  let  x  e  E(X)  and  x^  6  X^  .  Then 
x  -  E^x^  =  E(x-xh)  +  (E-E^Jx^  .  Consequently, 


S(x,Eh(Xh)}  < (  1 1 E  1 1 


(1.22) 


£(  I!  E|| 


!|x-xh||  +  ||E-Eh||x  i)xh||  ) 

h 

+  ||E-Eh||  x  )  It  x-xh||  +  ||  E-Eh||  x  ||x|| 


The  result  (1.3)  follows  by  taking  the  infimum  over  x^  e  x^ 
and  using  (1.1).  Q.E.D. 

It  remains  to  prove  Theorem  1.3.  We  first  quote  without  proof 
the  following  simple  result: 


Lemma  1.5.  Let  Y  and  Z  be  two  subspaces  of  X  with  the 
same  finite  dimension  and  let  P:  Y  -*■  Z  be  a  linear  operator 
such  that 

(1.23)  l|py-y||  1  \  lj y II  ,  y  e  y  . 

Then  P  is  bijective  and 

||P'1z||  <  2  ||z||  , 


(1.24) 


z  e  z  . 


Lemma  1.6. 

a)  If 

b)  If 


Let  Y  and  Z  be  two  subspaces  of 
5  (Y,Z)  <  1  ,  then  dim  Y  =  dim  Z  . 
dim  Y  =  dim  Z  <  » 

5 ( Y , Z )  <  5 (Z , Y)  [ 1-6 (Z,Y)  l"1 


X  . 


Proof.  For  a)  ,  See  Kato  [6,  p.  200]  .  For  b  ,  see  Kato  [7].  Q.E.D. 
In  this  rest  of  this  section,  we  suppose  that  the  hypotheses  of 
Theorem  1.3  are  satisfied.  Also,  c  and  hQ  will  denote  two 
generic  positive  constants  which  depend  on  r 


L-3S™jL-L.7.:.  5  (Al  E  (x)  ,AhlE.  (x.  ))  -  c6  (AI  E  (x)  ,Ah}  '  h  -  h0  * 

n  n 

Proof.  Let  x  e  E  (x)  ,  ||x|[  +  j|Ax||  =  1  ,  and  let  xh  e  Xh 

be  such  that 

||x-xh|j  +  ||Ax-Ahxh||  <  6(A|E(X)  ,Ah)  . 

We  have 

Hx"EhxtJI  +  11  Ax”AhEhxh^ 

■  il  2¥I  /.  (Rz  <a)  x-r2  <Ah) xh^ dz 

(1.25)  , 

+  il  2^1  /,  (ZRZ(A)X  -  ZRz(Ah)xh)dz|| 

<  (2*)-1  max  ||R  (A)x-R  (A.  )x.  II  /  (l+|z|)|dz|  . 

zer  2  z  n  n  r 

In  order  to  conclude  the  proof  of  the  lemma,  it  suffices  to 
estimate  ||  Rz  (A)  x-Rz  (Ah)  xh||  for  zer.  We  set  w  =  Rz(A)x  e  E(X) 


and  let  wh  e  Xh  be  such  that 
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(1.26)  ||w-wh||  +  ||Aw-Ahwhjj  <  6(A|E(x),Ah)  (  Ifw)!  +  |(Aw|(j. 
However , 

l|w||  +  ||Aw||  =  ||R2(A)x|1  +  ||  Rz  (A)  Ax|| 

1  ilVA)II  (  IN!  +  l!Axtl  ) 

*  II  r_  (a)  ||  <  c,  2  e  r  . 

Thus , 

(1.27)  ||w-wh||  +  ||Aw-Ahwh|j  £  e6  (aIE(x)  ,Ah)  * 

By  using  Theorem  1.1,  i.e.,  j|Rz(A^)||  ^  h  £  hQ,  zer, 

we  have 

nvvv-whif  *  |iRz(vK’(z'Ah)wh^ 11 

(1.28)  £  c  {  (I  |{  +  ||  (z-A)w- (z-Ah)wh||  > 

-  c6  (AIe(x)  ,Ah}  ' 

Consequently, 

HR2(A,x  "  Rz(Ah)xhl*  -  llw-whlf  +  Mwh"Rz(Ah,xh(i 

(1.29) 

£  ^A!g(3{)  /Ah^  *  Q.E. 

Remark .  Lemma  1.7  is  still  valid  if  the  hypothesis  PI)  is 
replaced  by  the  uniform  boundedness  of  R^A^)  on  r 


r 


'WgXSTrTlS’ 


Proof  of  Theorem  1.3.  From  (1.2)  and  (1.3)  of  Theorem  1.2  it 

A 

follows  that  lim  6  (E.  (X .  )  ,E  (x  ) )  =  0  .  Consequently,  by  Lemma  1.6  a), 

h-0 

dim  Eh(X^)  =  dim  E  (X )  =  m  for  h  <  h^  .  This  proves  that  for 
h  <  ^o'  Ah  ^  E  (X  ):  Eh^Xh^  Eh^Xh^  has  exact^y  m  eigenvalues 
^1  h',**,^m  h  rePeated  according  to  multiplicity. 

By  Lemma  1.6  b) ,  we  have  that 


(1.30) 


5lE(X)  ,Eh(Xh)  J  <  c  min{i(E(X),Eh(Xh))  , 

3tEh(Xh}  'E(X))}  '  h  <  hQ  . 


It  follows  from  Lemma  1.7  that 


6  (e  (X)  ,  Eh(*h))  i  C<S(AIE(X)  »AhlE  ,x  j) 

(1.31)  h  h 

<  c6(a|e(x) ,Ah )  . 

The  result  (1.31)  and  (1.2)  of  Theorem  1.2  yield  (1.4)  when 
substituted  in  (1.30)  . 

We  have  by  (1.1)  of  Theorem  1.2  that 

(1.32)  6^Ah^Eh(Xh)  '  -  c6(Ah'A);  h  <  h0  * 

Hypotheses  PI)  and  Lemmas  1.6  and  1.7  allow  us  to  conclude  the 
validity  of  (1.5)  . 

Now  let  u1#...,um  be  a  basis  for  E(X)  with  + 

II  Au .  1|  *l,i*l,...,m.  We  then  choose  u.  u  e  E.  (x.  ) 

*  iin  n  ti 


so  that 
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[|Ui-ui>h||  ♦  llAu^u.  bll  <  «l*k(x).AhlEh(!th,)  ' 

(1.33) 

x  3  1, . . .  ,m  « 

Next,  we  define  Ah:  E(X)  Eh(Xh)  as  the  linear  operator 
such  that  h,  i  =*  1, . . .  ,m  .  It  follows  by  (1.5) 

that 

(1.34)  ||u-Ahu||  +  ||Au-AhAhu|[  <  CYh  ||u||  ueE(X)  . 

By  Pi)  lim  y.  =  0  ,  so  by  Lemma  1.5  and  (1.34)  A,  is 
h-0  n  n 

a  bijection  whose  inverse  A^  is  uniformly  bounded  for  h  <  hQ 
Let  A  3  A|E(  }  and  ^  =  A^A^:  E(X)  E(x)  •  For  u  e  E(x) 
we  have 

||  (A-Ah)  u  ||  3  ||  A*1  (AhA-AhAh)  u  || 

d*35)  <_  c{  ||  (Ah-I)Aui|  +  !| Au-A^A^u ||  } 

<  CYh  ||u||  ,  h  <  hQ  . 

Consequently , 

(1*36)  l|A-Ah||  E(x)  £  CYh  . 

^  A 

Now  A^  has  eigenvalues  X^  h"**'^ra  h  an<*  A  **as  t*ie 
eigenvalue  X  of  algebraic  multiplicity  m  .  Also,  a  is  the 

A 

ascent  of  (X-A)  .  We  have  reduced  our  problem  to  the  matrix 
case.  The  results  (1.6)  and  (1.7)  of  Theorem  1.3  now  follow 
from  the  classical  perturbation  theory  for  finite  dimensional 
operators  (17,  p.  80-81].  Q.E.D. 
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2.  The  Approximation  of  the  Spectrum  of  a  Self adjoint  Operator. 

In  this  section,  we  suppose  that  x  is  a  Hilbert  space, 

A:  &  (K)  C  X  -*•  X  is  selfadjoint,  and  A.  :  X.  ♦  X.  is 

n  n  h 

selfadjoint  in  Xh  for  all  h  .  If  I  CJR  is  an  interval 
(finite  or  infinite)  ,  E ^ :  X  -*■  X  will  denote  the  spectral 
projector  of  A  relative  to  I  and  E,  _ :  X.  -<•  xu  will  denote 
the  spectral  projector  of  A^  relative  to  I  .  We  can  prove  the 
following  theorems: 


Theorem  2.1.  PI)  V  closed  intervals  I,J,  one  of  them 

bounded ,  I  n  J  =  ,  we  have 


lim 

h-0 


HEIEh,jl 


0 


Theorem  2.2.  Suppose  that  PI)  is  valid.  Let  J  Cl  where 
J  is  a  closed  bounded  interval  and  I  is  an  open  interval.  Then 

(2-D  lim  6(E  (X.),ET(X)J  =  0  . 

h-0  “  1 


Theorem  2.3.  Suppose  that  PI)  and  P2)  are  valid. 

a)  Let  I  be  an  open  interval  and  x  e  E^ (X)  .  Then 

lim<5((x,Eh>I<Xh))  *  0  . 

b)  If  X  e  a (A)  ,  then 


lim  dist  (x  ,a  (A.  ) )  ■  0 
h-*-Q  h 


/ 
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Remark  .  Theorem  2.3b  states  that  the  approximation  of  a (A) 
by  a (A^)  is  lower  semicontinuous . 

Theorems  2.1,  2.2,  and  2.3  have  been  proved  in  Section  3 
of  Descloux,  Nassif,  Rappaz  [4]  when  A  is  bounded  (recall 
our  remarks  after  Lemma  1.1).  Our  proof  ’-'ill  reduce  the 
unbounded  case  to  the  bounded  case.  We  shall  restrict  ourself 
to  the  proof  of  Theorem  2.2.  Theorems  2.1  and  2.3  can  be  obtained 
by  similar  arguments.  Furthermore,  for  the  sake  of  simplicity 
and  without  much  loss  of  generality  we  shall  suppose  that 
p  (A)  n  $  • 

Proof  of  Theorem  2.2.  As  mentioned  above,  we  suppose  there 
exists  a  6  ]R  with  a  S  p  (A)  .  We  introduce  the  function 
4>  (X)  =  (a-X)_1  and  set  B  =  <J>  (A)  =  R„  (A)  .  By  Theorem  1.1, 

*  <j>  (A^)  =  Ra(Ah)  is  well-defined  for  h  sufficiently 
small  and 

lim  ||  B-B  ||  =  0  . 

h-0  n  Xh 

For  an  interval  M  ,  FM:  X  -*■  X  and  F.  X,  ■*  X.  will  denote 

m  n,M  h  n 

respectively  the  spectral  projectors  of  B  and  Bh  relative  to 
M  .  As  a  first  case,  suppose  a  0  I  and  set  K  *  <t>  (J)  and 
L  =*  <J>  < I )  .  Then  K  and  L  are  respectively  a  compact  and  an 
open  interval  with  K  CL  .  From  Theorem  4  of  Descloux,  Nassif, 
Rappaz  [4 ] ,  we  can  conclude  that 

lim  6  (F.  „(X.  )  , F  (X))  =  0  . 
h^Q  h'K  h  L 


(2.2) 
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However,  the  result  (2.1)  follows  from  (2.2)  since  K(X^)  = 

Eh  J(Xh}  and  FL(X)  =  si(x)  * 

The  case  a  e  I  can  be  reduced  to  the  preceeding  one  by 
noting  that  we  can  find  compact  intervals  Jl,  J2  and  open 
intervals  II,  12  with  the  following  properties  (for  h 
sufficiently  small) : 


a  0  II  U  12  , 


Jl  C II  ,  J2  C  12  , 


Ej(X)  =  ES1(X)  9  ei2(X)  ' 
Eh,J(Xh>  *  Eh,Jl(Xh>  ®  Eh,J2IXh> 


Q.E.D. 


r 
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3 .  Application  to  Scalar  Elliptic  Boundary  Value  Problems. 

We  shall  need  the  following  notation  to  discuss  the  application 

of  our  theory  to  the  approximation  of  the  spectral  properties  of 

scalar  elliptic  boundary  value  problems.  Denote  by  Q  an  open 

set  in  lRn  with  a  smooth  boundary,  3^  .  As  usual,  we  denote  by 
2 

L  (fl)  the  Hilbert  space  of  square  integrable,  complex-valued 
functions  with  inner  product  and  norm 

(3.1)  <F,G>  =  fQ  FGdx  ,  !|F||2  =<F,F>  . 

We  denote  by  Hr(Q)  the  space  of  complex-valued  functions  whose 

distribution  derivatives  of  order  less  than  or  equal  to  r,r  a 

2 

nonnegative  integer,  are  in  L  (ft )  with  norm 

(3.2)  ||F,|2  =  l  |t  D°F||  . 

|a|  <  r 

We  wish  to  consider  the  spectral  approximation  of  the 
operator 


n  n 

Lu  =  -  7  D.  (a. .D.u)  +  7  a.  D ■  u  +  au  ,  x  e  a  , 

i,3*l  1  13  3  i=l  1  1 

(3.3) 

n 

=  l  a.  . v . D.u  *  0  ,  x  e  9Q  , 

i, j=l  13  x  3 

where  we  assume  that  L  is  a  uniformly,  strongly  elliptic 
operator  with  real-valued  coefficients  in  C° (Q)  and 
v  *  (v^,...,v  )  is  the  unit  exterior  normal  to  3D  .  We 
associate  with  L  the  continuous,  sesquilinear  form  on 
H1  (J2 )  x  H1  (£2)  , 
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(3.4) 


n  n 

i  <ai1D,^D.<:>  +  1  <a.D 

i,j=l  ]  1  i«l 


+  <a$,i|/> 


We  may  assume  (by  replacing  L  by  L  +  u  ,  a  e  IR  )  that  for 
some  b  >  0 

(3.5)  Re  B  (4> .  0 )  >  b||<M|2  ,  4>  e  H1(Q)  . 

Let  ^ h  C  H^(H)  be  a  family  of  finite  dimensional 
subspaces  parametrized  by  h  ,  0  <  h  _<  1  ,  with  r  a  positive 
integer  and  c  a  positive  constant,  independent  of  h  ,  such 
that  for  1  <  s  <  r  +  1  ,  u  e  Hs  ((2)  ,  we  have 


(3.6)  inf  (  Hu-xl!  +  h  ||u~xl|  ,  )  £  chs  |(u||  • 

x  e  ih 

Many  finite  element  spaces  are  known  to  satisfy  (3.6)  [2]  . 

We  define  the  operator  L^:  .S  h  -*•  3^  by  the  relation 

(3.7)  B(Uh,W)  =  <LhUh,W>  ,  W  6  ih  . 

Note  that  Lh  is  well-defined  since  is  finite  dimensional 

We  shall  study  the  approximation  of  the  spectral  properties  of 
by  the  spectral  properties  of  . 

Define  the  space 

H»(n)  =  {u  e  h2((2)  i  $u  =  o  for  x  e  an}  . 

2  2 

We  consider  L  as  a  closed  operator  from  L  ((2)  to  L  (Q)  , 

2  2 

i.e.,  X  =  L  (Q)  ,  with  domain  D(L)  =  (12 )  .  Also,  in  the 


-23- 


notation  of  Section  1,  we  set  .  It  follows  from 

( 3  -  6 )  that  if  u  e  C*(Q)  ,  then 

<5  (u,  .$.)  =  inf  j{  xa— x  II  -*•  0  as  h  -*■  0  . 

xe4h 

OO  —  2 

Since  C  (ft)  is  dense  in  L  (ft)  ,  it  is  clear  that  property 
P2)  is  valid. 

We  now  turn  to  the  verification  of  PI) .  To  that  end, 
let  P°:  L2(ft)  -*•  ^  be  the  L2(ft)  projection, 

—v 

(3.8)  <P°  u-u,W>  =  0  ,  W  e  Jh  , 

and  let  P^:  (ft )  -*■  ^  be  the  H^tft)  projection  defined  by 

(3.9)  b(p£u-u,w)  =  o  ,  we  Jh  . 

It  follows  from  (3.5)  that  P^  is  well-defined. 

The  following  estimates  for  P°  and  P^  are  well-known 
[2].  There  exists  c  <  <»  such  that  for  0  £  s  _<  r  +  1  and 
u  6  HS(!!)  , 

(3.10)  ||p£u-u||  <  chs  (| u (I  s 

and  such  that  for  1  _<  s  _<  r  +  1  and  u  e  Hs(i))  , 

(3.11)  ||pju-u|i  +h||pju-u|!1  <  chs||u  |js  . 

We  now  verify  that 

5  (Lh,L)  0  . 


(3.12) 
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For  CJ,  6  J,  ,  let  u  e  D(L)  satisfv 
n  n 

(3.13)  Lu  =  LhUh  . 

It  follows  by  elliptic  regularity  that 

(3.14)  Hull  2  <  c  |)LhUh||  . 

Now  it  is  easily  verified  that  =  P^u  .  Hence,  by  (3.11) 
and  (3.14)  , 

(3.15)  ||u-UhH  =  ||u-pju||  <  ch2  |I  u||  2  <  ch2  ||Lh0h||  . 

Since  Lu  =  L^U^  ,  we  can  conclude  that  for  each  S  , 

there  exists  u  e  D (L)  such  that 

(3.16)  l!u-Uh||  +  |(Lu-LhUh||  <  ch2(  ||uhj|  +  ||LhUh(i  )  . 
Hence,  we  can  conclude  that 

6 (Lh,L)  <  ch2  . 

It  is  well-known  that  the  spectrum  of  L  consists  of 
isolated  eigenvalues  with  finite  dimensional  generalized 
eigenspaces.  Let  X  e  c(L)  and  let  E (X)  be  the  generalized 
eigenspace  corresponding  to  X  .  Thus,  dim  E(X)  <  ®  .  We  shall 
show  that 

(3.17)  6  <L|e(x) ,Lh)  <  chr+1  . 

The  result  (3.17)  will  imply  that  the  conclusions  of  Theorem  1.3 
are  valid  with  Yh  ±  chr+1  . 


r 
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4 .  Application  to  the  Linearized  Shallow  Water  Equations. 

We  turn  now  to  the  description  and  analysis  of  an  approximation 
procedure  for  the  spectral  properties  of  the  operator  ,  T  , 
associated  with  the  linearized  shallow  water  equations.  Recall  that 

T(;,u)  =  (-7 •  u,-7; -fu-mRu)  ,  x  e  ft  , 

(4.1)  u  •  n  =  0  ,  x  6  3(2  , 

/  ^  dx  =  0 

(2 

2 

where  Q  C  ]R  is  a  bounded,  connected  open  set  with  smooth 
boundary,  3 a  ,  u  =  (u^,u2)  »  R  is  the  linear  operator 

Ru  =  (-u2,u^)  ,  n  is  the  exterior  normal  to  3a  ,  and  f  >_  0 
and  oj  are  real  constants  representing  friction  and  Coriolis 
terms.  We  assume  that  3(2  has  a  finite  number  of  connected  com¬ 
ponents,  •  Also,  assume  that  the  sets  3(2^  are  smooth 

2 

arcs  in  3R 

We  shall  need  the  definitions  of  the  following  spaces 
of  scalar- valued  functions: 

L*  ((2)  =  {w  e  L2(a)  |  /  wdx  *  0}  , 

(4.2)  Hi(n)  =  H1(n)n  L2(a/  , 

H^(ft)  =  {w  e  (a)  |  w(x)  -  0  for  x  e  3 f2 Q 

and  there  exists  constants, 

(ci)i_l  ,  such  that  w(x)  =  c.^ 
for  x  6  3Qif  i  =  1 , . . . , s }  . 

For  k  1  ,  set 

<4-3)  hJ(j2)  -  Hj(n)  n  Hk(n)  , 
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2  1 

We  give  these  subspaces  of  L  (ft)  and  H  (ft)  the  inner 
product  and  norms  defined  in  (3.1)  and  (3.2). 

We  also  define  the  following  spaces  of  vector-valued 
functions : 

L  (ft)  =  (u  =  (u1,u2) J u1,u2  6  L2(ft)}  , 

Hk(ft)2  =  (u  =  (ulfu2)|u1#u2  e  Hk(ft)}  for  k  >  1. 

With  these  spaces  wo  associate  the  inner  product  and  norms 
<u,v>  =  <u1,v1>  +  Cu2,v2>  , 

II  u  ||2  =  <3,3  >  , 

for  u  =  (u^,u2),  v  *  ^vi'v2^  6  L  ^  ' 

IIs  Ilk  *  +  II“2II k 

for  3  =*  (u1,u2)  e  Kk  (ft )  2  . 

Finally,  for  functions  (?,3)  ,  (£,v)  e  L2(ft)  x  L2(ft)2 
we  define  the  inner  product 

<((?,U),  (CfV)>  =»  +  <U,V> 

and  norm 

|!  (?»u)  ||  2  a  <(c  ,u)  ,  (£,U)  > 

and  for  functions  (?,u)  6  H*(ft)  x  L2(ft)2  we  define  the  norm 

-  lUIlf  +  l|3||2 


l|(«.»>llx 


/ 
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(4.6) 


,RVipy  =  0. 


<t>  e  h*  (0) ,  i/>  e  (u)  . 


Hence,  (4.4)  is  an  orthogonal  direct  sum  decomposition. 

Let  (7*1^  >  0  <  h  <_  1  ,  be  a  family  of  finite  dimensional 
subspaces  of  H^(ft)  parametrized  by  h  .  Set 

(4.7)  =  ?>?h  O  H*(Q)  , 

(4.8)  7^*  =•  7>lh  r\  Hi  (ft)  . 

We  assume  that  there  exists  a  positive  constant  c  , 
independent  of  h  ,  and  a  positive  integer  r  such  that  for 
1  1  k  1  r  +  1  and  w  e  Hk(fl)  , 

(4.9)  inf  (  ||  w-x  ||  +  h  ||w-xl!  x  )  1  chk  ||w||k  , 

xe*; 

and  for  1  _<  k  ^  r  +  1  and  w  e  H^(!5)  , 

(4.10)  inf  e  (  II  w-x  ||  +h||w-x||1)  <  chk  ||w||  . 

X  en£ 


We  also  wish  to  assume  that  the  spaces  {7*^}  satisfy  the 
"  inverse  property"  that  there  exists  a  positive  constant 
c  ,  independent  of  h  ,  such  that 
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We  shall  use  the  spaces  of  vector-valued  functions 

71 h  =  {uju  =*  -v*  +  rw,  $  e  Tn^  >  *  e  7>ip 

It  follows  from  (4.5)  ,  (4.9),  and  (4.10)  that  for 
0  £  k  <  r,  u  e  Hk(Q)2  , 

(4.11)  inf  j|u-x!i  <  HujJ  .  . 

X^h 
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•  c 
We  shall  need  the  result  that  if  f  0  77?  h  ,  then 

(4.12)  7?  e  ^?h  . 

It  is  easy  to  see  that  (4.12)  is  valid  if  contains  the 

constant  functions. 

In  order  to  define  our  approximation  procedure,  we 
introduce  the  continuous ,  sesquilinear  form  B ( - , • )  on  X  *  X 
[H*  (ft)  x  L2(ft)2]  x  [«£(&)  x  L2(a)2]  by 

B((S,u),  (£,v))  =  <u,7£>  -  ,v> 

(4.13)  -  <fu+wRu,v  >  , 

« ,u) , (c,v)  e  h* (n)  x  l2 (Q) 2  . 


We  regard  T  in  (4.1)  as  an  unbounded,  closed  operator, 
T:  H*(S1)  x  l2(G)2  *  H^(n)  x  l2(G)2  ,  with  domain 

d (t)  =  { te  ,u)  |c  e  hJ («) 

(4.14)  u  e  L2<a)2  ,7  •  u  e  h1^)  , 

u  •  n  *  0  on  30  }  . 


Thus,  if  te  , u)  e  D (T)  , 


(4.15) 


Blfc  ,u)  ,  (C,v)j  *  <T(C,u),(£,v)>  , 


(£,v)  e  Hj(fl)  x  l2(Q)2 


We  now 


define  Th:  7*1  h  *  ♦  7*7^  x  bY 


(4.16) 


B  [Z  ,U)  ,  (Y  ,V) )  «  <Th(Z,?J)  ,  (Y,$)  >  , 

(Y,$)  e  7fi*  x  7^ 
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The  operator  T^  is  well-defined  since  the  spaces  x 

are  finite  dimensional. 

We  now  show  that  the  spectral  properties  of 

approximate  those  of  T  by  verifying  PI)  and  P2)  in  this 

case  for  X  *  a*(n)  *  L2(ft)2  and  X.  =  Tti*  *  W  .  The 
*  h  h  h 

verification  of  P2)  follows  directly  from  (4.9)  and  (4.11). 

To  see  this,  we  note  that  for  (S,u)  e  L2  (ft)  O  C* (Q)  *  c'”(?2)2 
it  follows  from  (4.9)  and  (4.11)  that 

$((C,u)  ,Xh]  ■*  0  as  h  -  0  . 
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9  _  _  30  —* 

The  validity  of  P2)  now  follows  from  the  density  of  L*  (71)  0  c  (.1)  * 

C*  (Q )  2  in  H*  <fl)  x  L 2  (ft)  2  . 

We  next  turn  to  the  verification  of  PI)  .  In  order  to 
prove  this  result,  we  introduce  and  analyze  some  projection  opera¬ 
tors.  First,  we  define  the  H*  (12 )  projection  p*-.  ni<a)  - 
by 

(4.17)  <7 (P*w-w) ,7X>  -  0  ,  x  e  ^  , 

and  the  (S3 )  projection  P°:  H^(0)  -*•#!■£  by 

(4.18)  <7  (P*w-w)  ,7x>  *  0  ,  x  e  '»l\  • 

It  is  well-known  [2]  that  there  exists  c  <  <*>  such  that  for 
1  <  k  <  r  +  j.  and  w  6  H*(!))  , 

(4.19)  ||P*w-w||  1  <  ch*'1  ||w|lk 

and  for  1  <  k  <  r  +  1  and  w  e  Hk (ft )  , 

—  —  c 

(4.20)  llP^w-wlli  <  chk_1  i! w jj R  . 

It  follows  from  (4.11)  that  if  <$.  :  L2(Q)2  ■*  71  is  the 

n  a 

2  2 

L  (ft)  projection  defined  by 

(4.21)  <  $hu-u,x>  *  0  ,  X  e  7lh  , 

then  there  exists  a  constant  c  <  »  such  that  0  <  k  <  r 

V  O 

and  u  e  H  (Q)  implies  that 

(4.22)  ll<5hu-u||  <_  chk  (i  ui|  k  . 


It  is  easily  checked  that 
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(4.23)  Q.u  =  -7P*Su  +  RVP^Hu  . 

ha  h 

Finally,  define  the  L2(^)  projection  Q^:  L2(H)  -  7*1  *  by 

(4.24)  <Qhw-w,x>  =0,  X  6  . 

It  follows  from  (4.9)  and  the  inverse  property  of  that  there 

}c 

exists  c  <  ®  such  that  if  l£k£r+l  and  w  6  H*(ft)  ,  then 

(4.25)  H^hw_w!ll—  Chk  ^  ' 

The  following  lemmas  will  help  us  analyze  . 

Lemma  4.1.  a)  Let  u  e  L 2  (S2 )  2  and  x  6  .  Then 

(4.26)  <u-$hS,7X>  =  0  . 

b)  Let  5  e  ai(Q)  and  X  6  .  Then 

(4.27)  <7(5“PgO  »X>  -  0  • 

Proof.  The  proof  follows  easily  from  (4.6)  and  the  definitions 
of  P*  and  <5^  .  Q.E.D . 

Now  let  <?h  *  P*  ®  3h:  1*1(0)  x  L2(ft)2  -  7)?*  x  . 

Lemma  4.2.  If  (?, u)  e  D(T)  and  (Y,V)  6  7$*  x  ^  ,  then 

B((c,u)  -  ^h(C,u),(Y,$))  - 

(4.28)  <T(c,u)  -  Th^h(C,u)  ,  (Y,^?)  > 

*  -  <^a>R(u-$hu)  , V  > 

Proof.  The  proof  follow  directly  from  the  definition  of  B  and 
Lemma  4.1.  Q . S . 0 • 
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Hence,  it  follows  from  Lemma  4.2  that  if  Q.  =  Q.  ®  Q.  , 

h  h  h 

(ou)  e  D(T)  ,  and  (Y,V)  e  x  ^  ,  then 

(4.29)  <QhTU/U)  -  Th^h(;,u)  ,  (Y,V)> 

=  -  <wR(u-$hu)  ,V>  . 

Theorem  4.1.  <S(T,_,T)  <  ch. 

-  n  — 

Proof.  Let  (Z,U)  e  7^*  x  and  denote  Th(Z,U)  by  (D,E)  . 

We  must  fincf  (S,u)  e  D  (T)  such  that 

(4.30)  ||  (Z,U)-(C,u)  ||x  +  ||Th(Z,U)-T(?,u)  1^ 

<  Ch(  ||  (Z  ,U)  |^  +  ||Th(z,u)  |y  . 

We  choose  (C,u)  as  follows.  Let  $  e  H*(ft)  be  the 
solution  to 


(4.31) 


A<f>  =  D  * 


x  6  n  , 


=  0  x  6  3ft 

®n 


3 

where  ^  is  the  exterior  normal  derivative.  (Note  that 

/  D  dx  *  0  since  D  €  L^* (ft)  ) .  Now  set  ?  =  Z  ,  u  =  — + 

ft  _  1+1 

R7HU  .  It  is  clear  that  since  Z  e  H*(ft)  ,  7  •  u  =  — A<J>  *  -D  6  H  (ft)  , 

u  •  n  »  -  “  0  for  x  e  9ft  that  (?,u)  e  D  (T)  . 

It  follows  by  the  definition  of  $  that  P*$  *  S(5  .  By  elliptic 

regularity  for  (4.31), 

IU  II  2  1  c  II D  ll  1  c  H D  ||  x  . 


Hence,  we  can  conclude  from  (4.19)  that 
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(4.32) 


||V(Su-S0)li  =  ||7(«-P*0)|J 

<  ch  ||$||  2  1  ch  !! D I!  x  . 


So, 

II  (Z,U)  -  U,u)  II  x 

(4.33)  =  j|  V(Su-SU)  l!  £  ch  ||d|| 

<  ch  ||Th(Z,U)  Hx  . 

» 

Denote  T(S,u)  by  (d,e)  .  It  follows  from  (4.31)  that 

d  =  D  .  Since  @ .  U  ,u)  =  (2,U)  ,  we  obtain  from  (4.29) 

n 

(4.34)  <Qhe-E,V>=  <«R(7*-7P*<M  ,V>  . 


Thus,  if  we  set  V  =  Qhe  -  E  in  (4.34)  and  use  the  Cauchy-Schwarz 
inequality  we  can  derive  the  estimate 


(4.35) 


|QhS-E||  <  c  ||74>-7P**!|  <  ch 1 1  1 1  2 

<  ch  ||  D  ||x  <  ch  ||Th(Z,U)  Hj 


We  must  now  estimate 


(4.36) 


e-Q,  e||  *  _^inf  ||e-xl| 

xe77. 


By  the  definition  of  T  , 

e  =  -7  C  -  fu  -  wRu 
(4.37)  *  -  7Z  -  f  (-V4>  +  RVHU) 

-  u){-R74>  -  VHU)  . 


However , 
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(4.38) 


-7Z  -  f  R7H$  +  ^>7  HU  e 


and 

(4.39)  i|  f7$  +  UR7$  ||  L  <_  c  ||  D||  x  . 


So, 


||e-$.e||  =  inf  ||  f7$  +-uRV<p-x  il 

x^l, 

(4.40)  1  oh  ||  £7*  +  uR7*  ||  , 

<  ch  ||  D !| x  _<  Oh  ||Th(Z,U)  ||x  . 

The  triangle  inequality,  (4.35)  and  (4.40)  yield 


(4.41) 


||T(C,u)  -  Th(Z,U)||x 


*  l|e-E||  <  ch  ||Th(Z,U)  ||x  . 

Thus,  (4.30)  follows  from  (4.33)  and  (4.41).  Q.E.D. 

Now  let  X  e  a(T)  be  an  isolated  eigenvalue  of  T  with 

finite  algebraic  multiplicity  m  .  Let  E ( X)  be  the  generalized 

eigenspace  corresponding  to  X  .  It  follows  from  Theorem  1.3 

that  we  are  interested  in  an  estimate  for  5  fT  I  ,  ,T.  }  . 

v  E  (X)  h' 


Theorem  4.2.  Let  X  6  a(T)  be  an  isolated  eigenvalue  of 
T  and  let  E(x)  be  its  associated  generalized  eigenspace. 
Assume  dim  E(X)  =  m  <  »  ,  and  assume  that  U  ,u)  S  E(x) 
implies  that  ?  e  H*+1(ft)  ,  u  e  Hr(n)2  ,  7*u  e  Hr+1(n). 
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Then  there  exists  c  <  ®  ,  independent  of  h  ,  such  that 

(4.42)  6(tIe(X)  ,Th^  -  Ch  ' 

Proof.  Let  (;  ,u)  6  E(X)  ;  Set  (2,0)  =  <^hU-u)  •  It 

follows  from  (4.19)  and  \ 4 . 22)  that 

(4.43)  ||  (C,u)  -  rfhU,u)  ||x  £  chr(  II  :||r+1  +  • 

Also,  setting  (Y,V)  =  QhT(;,u)  -  V*h(5'S)  in  (4.29)  and 
using  the  Cauchy-Schwarz  inequality  we  obtain 

||QhT(;,u)  -  Th  ^h(;,u)llx  =  !lQhT(;fu) -Th^h(;,u)  I! 
(4‘44)  <  c  ||u-Qhu|j  <_  chr  ||u|1  r  • 

Finally  from  (4.22)  and  (4.24)  we  obtain 

(4.45)  ||T(5,u)  -  QhT(5,u)  ||x 

£  <=hr(  IUIIr+1+  ||u||r+  !|V-ui|r+1) 

The  result  of  the  theorem  now  follows  from  (4.43),  (4.44), 
(4.45)  and  the  finite  dimensionality  of  E (X)  .  Q.E.D. 

Finally,  we  note  that  our  results  can  be  combined  with 
the  results  in  [3]  to  obtain  optimal  order  estimates  on  the 
convergence  of  eigenvalues .  Assume  that  the  conditions  of 
Theorem  4.2  are  valid.  Let  E(X*)  denote  the  generalized 
eigenspace  associated  with  the  eigenvalue  X  of  T*  .  We  note 
that  it  is  well  known  that  dim  E  (X*)  =  m  and  that  X  is  an 
isolated  eigenvalue. 
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Let 

wh  =  6 (E (X) , E (Xh) )  , 
w*  =  6 (E (X* ) ,E(xJ>)  . 

Let  a  be  the  ascent  of  X  .  Then  it  follows  from  Proposition 
3.2  of  [3]  that  for  h^  sufficiently  small 


max  |X-X^  h|a  <_  c  w^w*  , 


i — 1 / .  •  •  ,m 
m 


X  -  —  j  X .  .  I  <  c  w,  w*  ,  h  <  h_ 
m  i,h'  —  h  h  —  0 


Under  the  hypotheses  of  Theorem  4.2  we  have  proven  that 
wh  _<  chr  .  Now  assume,  in  addition,  that  U,u)  e  E(x*) 
implies  that  c  e  H*+1(n)  ,  u  e  Hr ( ^ ) 2  ,  ?  *u  6  Hr+1(ft). 
Then  we  can  conclude  from  applying  the  above  arguments 
to  the  adjoint  problem  that  w*  £  chr  .  Hence,  under  the 
above  conditions  we  obtain  the  optimal  order  eigenvalue 
estimates 


max 

i*l» . . . ,m 


'  J'xi,hl  i  ch 


2r 


m 


l 

i*l 


<  ch 


2r 


h  <  h 


0  * 
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